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A discrete time crystal (DTC) is a robust phase of driven systems that breaks the discrete time 
translation symmetry of the driving Hamiltonian. Recent experiments have observed DTC signa¬ 
tures in two distinct systems. Here we show nuclear magnetic resonance (NMR) observations of 
DTC signatures in a third, strikingly different system: an ordered spatial crystal. We use a novel 
DTC echo experiment to probe the coherence of the driven system. Finally, we show that inter¬ 
actions during the pulse of the DTC sequence contribute to the decay of the signal, complicating 
attempts to measure the intrinsic lifetime of the DTC. 


Periodic driving of a many-body system may lead to in¬ 
teresting out-of-equilibrium states of matter. A recently 
described phenomenon that has received much attention 
is the discrete time-crystalline phase, in which a driven 
system spontaneously breaks the discrete time transla¬ 
tion symmetry of its underlying Hamiltonian (also known 
as a Floquet time crystal, or a 7r-spin glass) [IHBj . Follow¬ 
ing a proposal for how to realize a discrete time crystal 
(DTC) [7], two experiments in quick succession showed 
evidence of DTC order. The first used trapped ions 
[SJ, and closely followed original theoretical models for 
DTC behavior, specifically working in a regime thought 
to be conducive to many-body localization (MBL) [9]. 
The second used diamond nitrogen vacancy (NV) cen¬ 
ters, which introduced several novel features: many spins, 
three-dimensional geometry, and long-range dipolar in¬ 
teractions m- Although systems with these traits are 
considered unlikely candidates for MBL, sources of dis¬ 
order did exist in that experiment, leaving some uncer¬ 
tainty as to the role of MBL in the observed signatures 
of DTC order [TTUTi] . 

In this Letter, we report the observation of discrete 
time-crystalline signatures in a system expected to be 
even further from the MBL limit than prior systems: an 
ordered spatial crystal. We observe robust oscillations at 
half the drive frequency (“DTC oscillations” for brevity) 
across orders of magnitude in interaction time (Figs, [l] 
and [2j . We also study the decay mechanism of the DTC 
oscillations, with two results. First, we show by generat¬ 
ing a time-reversed DTC echo that the density matrix is 
more coherent than the original DTC sequence reveals. 
Second, we show that the effect of interactions during 
the nonzero pulse duration of the DTC sequence limits 
our ability to observe the intrinsic lifetime of the DTC 
oscillations. 

We study the 100% occupied crystal lattice of spin-1/2 
31 P nuclei in ammonium dihydrogen phosphate (ADP), 
with chemical formula NH4H2PO4. ADP is an ionic, 
tetragonal crystal that also contains spin-1/2 4 H nuclei 
(99.98% abundant) and spin-1 14 N nuclei (99.64% abun¬ 
dant) [T51 . 

In the rotating frame pa. the secular internal spin 


( a ) 20 



-20 1 . 

0 20 40 60 80 100 120 

Floquet cycle (N) 



0.4 0.5 0.6 


v 


FIG. 1. After preparing the 31 P spins in a +z state, we 
measure the signal after the pulse sequence {r—X n } N . Green 
dots are for N odd, blue are for N even, and lines are to guide 
the eye. (a) For r = 12.5 ps and 0 = 0.9957T, the response is 
an oscillation at half the drive frequency, corresponding to a 
single peak in the Fourier transform at normalized frequency 
v = 1/2. (b) At the same r = 12.5 ps but with 9 = 1.0547T, 
the signal exhibits a beat pattern, corresponding to a splitting 
of the Fourier peak, (c) For increased r = 392.5 ps and 
6 = 1.060tt, the sharp peak at v = 1/2 is restored, despite the 
deviation of 9 from 7r. 


Hamiltonian for 31 P is 

nj _ o/P _j_ o/P,P 1 o/P,H . m_/P,N /-I \ 

ri int — riz "T Tt Z z + Tt zz • UJ 

Here, ~ S2 / Zt is a single-spin Zeeman term due to a 
uniform resonance offset El: and fizz describes the var¬ 
ious secular dipolar couplings between spins (where we 
omit terms that do not involve 31 P) [19j. These dipo¬ 
lar coupling terms are summarized in Table [TJ with typi- 

























2 


TABLE I. Terms in the Hamiltonian involving 31 P nuclear 
spins in ADP. The dipolar coupling constant for a pair of 

spins is Bij = ^ \ [1 — 3cos * 2 ( dij )], where Oij is the angle 

between the internuclear vector r)j and the z axis (defined 
by the static external field), yo is the vacuum permeability, 
and 7 i and 7 , are the nuclear gyromagnetic ratios of the two 
spins. {1$, S/>, R 4 ,}, (j) = x, y,z, are the spin operators for 
{ 31 P, 1 H, 14 N}. I 11 the presence of the external Ho = 4T field, 
the Larmor frequencies cuo/27t = 'yHo/2iv of 14 N, 31 P, and 4 H 
nuclei are 12, 69, and 170 MHz respectively. 


Hamiltonian term Typical interaction freq. (Hz) 

n*i p = -h-Ij) 11 /P,p /27t = 508 

H F i s = /l'j(2/,„.SV v ) W F ’ H /2n = 3500 

^ = Eij Bfi (2/* Rz, ) W p ’*/2n = 97 


cal interaction frequencies W/2n from the numerics that 
best match our experiments m ■ Note that the dipo¬ 
lar coupling is a long-range interaction (oc r~ 3 in a 
three-dimensional system). During a strong rf pulse of 
phase 4>, an external pulse Hamiltonian H r f = 
is added to Eq. ([!]). Since wi/27t ~ 68 kHz is so large, 
the internal Hamiltonian terms are often ignored during 
the pulse action (the delta-function pulse approximation) 

Q3H2]. We revisit this approximation later. 

There are four primary differences between the spin 
Hamiltonian in this system and prior models. First, the 
spin Hamiltonian is exceptionally ordered. The Zeeman 
term is extremely uniform in space, and the symmetry 
of the ADP crystal causes each 31 P nucleus to experi¬ 
ence an identical set of dipolar couplings from other 31 P, 

4 H, and 14 N jTBj. Second, rather than having only an 
Ising-type term I z , I Zj , the dipolar coupling between the 
31 P nuclei includes “flip-flop” terms I Xi I Xj + lyily = 
(itIJ + ir^t)/ 2- We exploit this difference to make 
echoes, as we describe below. Third, while the detected 
spin-1/2 31 P spins interact among themselves, they also 
interact with two other spin species that are not directly 
affected by the repeated pulses: the spin-1 14 N (with a 
weaker coupling) and the spin-1/2 4 H (with a stronger 
coupling). Fourth, we are able to selectively use high 
power continuous wave (cw) decoupling to “effectively 
remove” the 4 H from the spin Hamiltonian; we label this 
111 H off.” Experiments where this technique is not used 
are labeled ul H on.” 

To do these experiments, we first prepare the 31 P spins 
in a weakly polarized +z state (not a pure state) by 
cross polarizing the 31 P spins with the bath of room- 
temperature 4 H spins, producing an initial +y polariza¬ 
tion of 31 P, then applying a final pulse to rotate the mag¬ 
netization into z [16] . While this procedure improves the 
initial polarization by a modest factor 7 h/7p ~ 2.5, it 
significantly increases our ability to quickly repeat ex¬ 
periments, as the relaxation time T-/ 1 = 0.6 s of the 4 H is 
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FIG. 2. To explore the onset and robustness of oscillations at 
half the drive frequency, we examine the crystalline fraction 
/ using decoupling to turn the 4 H off (left) and on (right), 
[(a) and (d)] The crystalline fraction varies smoothly as a 
function of 9/n, fitting to Gaussians. Crystalline fractions 
from r = 12.5 (open triangles) and 392.5 ys (closed circles) 
are shown. The horizontal dotted line shows / = 0.1. [(b) and 
(c)] Cutoff at / = 0.1 from the Gaussian fits, corresponding 
to the region in 9 and r within which we observe persistent 
DTC oscillations. With the 4 H turned off, the boundary of 
this region shows structure near r = 1 ms (see [16] for further 
discussion). We also show \9 — 7 r| = W F ' F r [black dashes, (b) 
and (c), and (e) and (f)], where 1 T P ' P is the typical interaction 
scale of the 31 p_ 31 p coupling (Table |l}. For r > 10 ms, the 
decoupling power begins to heat the tank circuit, skewing 
results, and preventing exploration of the r > 10 ms region 
for 4 H off. [(d)-(f)] Without decoupling, we can avoid heating 
and explore further in r with 4 H on. I 11 (f), data span the 
range 0.03 < lT p p r < 3200 radians. 


much faster than that of the phosphorus, Tf* = 103 s. 

Beginning with the enhanced 31 P spin polarization, we 
apply the pulse sequence {r — Xg} N , where r is a wait 
time, and the pulse Xg is a rotation about +x by angle 
6. After N such repeated Floquet cycles, we apply a final 
X ^/2 to generate a nuclear magnetic resonance (NMR) 
signal that reveals the z magnetization after N cycles. 
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Unless otherwise stated, we apply 1 H decoupling [22] 
from the end of cross polarization through the readout. 
After reinitializing the starting state over time ss 5T| H . 
the pulse sequence is repeated for the next TV, from N = 1 
to N = 128, producing a discrete time signal S(t) with 
Fourier transform S(u). In practice, the internal Hamil¬ 
tonian [Eq. 0] acts not only during wait time r, but 
also during the pulse, which is of duration t p ~ 7.5 /zs. 
The total Floquet period is T = r + t p , with drive fre¬ 
quency vp = 1 /T; below, we examine the spectrum S(u) 
as a function of the normalized frequency v = v/vp. 

At small r and 9 ss w, the system responds with trivial 
oscillations at half the drive frequency, i.e., at v = 1/2 
[Fig. 0 a)]. As 9 is adjusted away from n, the response 
beats with frequency determined by 9 — tt [Fig. 0b)], 
also an expected result. However, with 9 ^ tt and in¬ 
creased values of r, the system returns to robust v = 1/2 
oscillations, a signature of DTC behavior 0 [Fig- 0c)]. 

We explore this behavior as a function of r and 9 
by examining the “crystalline fraction” / = \S(u = 
1/2)| 2 /|S'(i>)| 2 , as defined in [ID] , except that we 
Fourier transform our entire S(t) rather than just a late¬ 
time window (see [ID] for further discussion of this point). 
For given r, we fit the crystalline fraction as a function 
of 9 to Gaussians with good results [Fig. 0a)]. Note 
that the precise shape of f(9) does depend on the choice 
of Fourier transform window size, as described in m- 
These Gaussians reveal a region within which the robust 
v = 1/2 oscillations are detectable (the “DTC region”), 
and outside of which there are diminished or split Fourier 
peaks. We may visualize the DTC region [shown in Figs. 
0 b) and0c)] by choosing a cutoff / = 0.1. The width 
of the DTC region increases with r at short r, fluctuates 
slightly near r = 1 ms (discussed further in mi). and 
then approaches a steady value at long r. 

Although we explore multiple decades with 1 H off, 
heating of the NMR tank circuit (from cw decoupling) 
imposes an experimental barrier for r > 10ms [IB]. In 
order to explore further in r, we turn off the decoupling, 
allowing the X H to act. We prepare the initial state in the 
same way as described above, but then apply pulses to 
the 31 P only [see Figs. 0d)0f)]. In this 1 H-on case, the 
DTC region approaches its maximum width faster as a 
function of t, and shows little to no fluctuation in width 
over many decades in r. The width appears to decrease 
slightly when the experiment time approaches Tj p , the 
relaxation time of the 31 P spins. 

While these robust oscillations for many r are distinc¬ 
tive signatures of a discrete-time-crystalline phase, the 
lifetime of the state is of particular interest. To explore 
reasons for the decay of oscillations during the DTC ex¬ 
periment, we begin by considering a trivial decay mech¬ 
anism for noninteracting spins m- If we start with a 
net magnetization vector aligned along z, then an X„+ e 
pulse leaves a component of the magnetization along the 
x axis proportional to cos(e). If the component of the 



FIG. 3. A decay envelope cos N (e) (solid black line) bounds 
the magnitude of our DTC oscillations (closed red triangles) 
for long r and large e. To show that the observed decay 
is not dominated by an irreversible process, we devise an ap¬ 
proximate “unwrapping” sequence to create an echo above the 
classical decay envelope. After N cycles of {r — Xg} (closed 
red triangles), we switch to the second part of the echo se¬ 
quence (see text) and increment N'. The results of this N' 
sequence are shown for r = 192.5 p,s, where clear echoes are 
observable for both 9 = 1.087T (a) and 9 = 1.167T (b), with 
N = {2, 6, 10} (open blue circles, green triangles, and yellow 
diamonds). The corresponding arrow and filled marker show 
the predicted time of the DTC echo maximum, i.e., where 
N' = N. For the case of the DTC sequence at 9 = p (blue 
dots), we observe decay where none is expected in the perfect, 
delta-function pulse model; the echoes do not rise above this 
9 = tv decay envelope. We examine this effect in Fig. [I] 


magnetization in the transverse plane after the pulse is 
assumed to dephase and be lost during the subsequent 
r due to some variation in local fields, we expect each 
Floquet cycle to reduce the observable magnetization by 
cos(e), leading to an exponential decay cos N (e) after N 
Floquet cycles. In practice, our signals seem to stay at or 
below the barrier imposed by this decay rate (compare 
Fig. 0 closed red triangles versus black solid line). 

To test whether magnetization is irreversibly lost dur¬ 
ing the decay of the DTC oscillations, we devise a “DTC 
echo” sequence that attempts to “undo” the forward evo¬ 
lution of the system during the original DTC sequence. 
In the original sequence {r — Xg} N , the system evolves 
under TAnt for time r, followed by a pulse rotation. To 
undo the rotation, we simply apply a pulse rotation in 
the opposite sense, Xg. To undo the effect of "Hint dur¬ 
ing r, we first limit our attention to the homonuclear 
interaction , and we borrow a technique from the 
so-called “magic echo” experiment in NMR [191 23] , ex¬ 
ploiting the fact that a strong rf pulse along y reduces 
H F ' F to an effective internal Hamiltonian term — h'H F,F , 

z yy 

where we define = Y,^>i B fj( ~~ £ ’ Ij) PH- 
By allowing — \ B yy > to act for time 2r and rotating the 
state with surrounding X n / 2 pulses, the net phase evo¬ 
lution becomes opposite the original forward 

dipolar evolution The final DTC echo sequence 
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is 

{r - X e } N - (A w/2 - {X e - Yz} n ' - X w/2 ), (2) 

where $ = u>\2t for a strong rf pulse of duration 2 r. 
In the language of the more conventional Hahn spin echo 
sequence, the first part (N blocks) of this sequence gen¬ 
erates the “free induction decay,” while the second part 
(the rotated N' blocks) generates the signal “after the n 
pulse” [T32 [24]. If the second part of this sequence can 
undo the first part, then we expect an echo to appear 
when N' = N. In Fig. [3j we show prominent echoes in 
the magnetization for various values of N. The detection 
of these DTC echoes is evidence that the original DTC 
sequence drives coherence to normally unobservable parts 
of the density matrix, without being irreversibly lost due 
to interactions with an external thermal bath. 

Note that even at 9 = ir, the signal decays (Fig. [3]), 
which should not happen in the delta-function pulse ap¬ 
proximation. Moreover, this 9 = ir decay envelope ap¬ 
pears to limit the recovery of signal for 0 ^ ir, where we 
see that the DTC echoes are unable to cross this bar¬ 
rier (Fig. [3]). This may be also responsible for the echo 
peaks occurring earlier than their expected locations. An 
obvious mechanism that could cause this is pulse imper¬ 
fections; however, our studies have shown their effects to 
be too small to explain the observed decay [16]. 

In order to account for this decay, we reexamine the 
previously ignored effects of 7Ant during the nonzero du¬ 
ration pulses. Experimentally, we explore this by us¬ 
ing modified versions of the DTC sequence, using differ¬ 
ent sets of pulse phases at 9 = tt. Defining {a,/3} = 
{r — ol^ — t — /3tt} n , we expect little or no difference 
between {A, A}, {A, A}, and {X, Y} for ideal delta- 
function pulses. However, Fig. |4][a) shows that while 
{A, A'} and {A, A} produce very similar results, {A, A} 
yields a dramatically extended lifetime. 

This result can be qualitatively understood when ac¬ 
counting for the nonzero pulse duration t p . During a 
pulse of phase A (A), an effective internal Hamilto¬ 
nian — §?^. p (— \'H- T y’y) acts f° r ti me tp H3 I2B- For 
{A, A'}, a term proportional to is present in the 

average Hamiltonian; this term is for {A, A}. 

However, for {A, A} this term is instead proportional 
to — -I -Hy^) = +^T~L?f JSj, which commutes 

with the initial density matrix. We explore the effect of 
nonzero pulse duration further in |16l . 

The results we have shown are strikingly similar to 
those described by Zhang et al. [8] and Choi et al. CD], 
even though the spin Hamiltonian for our system is dif¬ 
ferent in interesting ways. We also explore a very large 
region in the (9,t) parameter space, and observe DTC 
signatures across a remarkably broad range in r. This 
occurs despite the conventional wisdom that our experi¬ 
mental conditions should be even less conducive to MBL 
than all prior DTC experiments jJJ 112) . DTC oscilla- 



Time / T* 


FIG. 4. Experimental pulses are of nonzero duration, which 
allows the internal Hamiltonian to act. These actions, while 
small, produce significant effects, especially after many Flo- 
quet cycles. At r = 20 /rs, pulse sequences {A, A'} (black open 
squares) and {A, A} (red open circles) decay much faster than 
the sequence {A, A} (green open triangles). Signals are sam¬ 
pled after each repeating block {a,/?}. Here, our T* = r + t p . 


tions like those observed in this experiment can some¬ 
times be explained as prethermal phenomena, but our 
system seems to violate conditions for prethermalization 
HQ US]. The DTC echo opens a new window into the 
physics of these driven systems, and helps to clarify the 
nature of the DTC oscillations. Further, we have shown 
that the interactions during the pulse of the DTC se¬ 
quence contribute to the decay of the oscillating signal. 
This is a practical barrier to measuring the intrinsic life¬ 
time of the DTC, which should be taken into account in 
future studies. 

In this work, we exploited both the long coherence 
times of our sample, and our ability to use NMR pulse 
sequences to edit the spin Hamiltonian. This suggests 
that NMR can be a useful probe of the physics of out-of- 
equilibrium, driven many-body systems. 

We thank C. W. von Keyserlingk, V. Khemani, C. 
Nayak, N. Yao, and M. Cheng for helpful discussions. 
We also thank C. Grant and D. Johnson for help in con¬ 
structing the NMR probe, K. Zilrn for recommending the 
ADP sample, and S. Elrington for assistance with imple¬ 
menting cross polarization. This material is based upon 
work supported by the National Science Foundation un¬ 
der Grant No. DMR-1610313. R.L.B. acknowledges sup¬ 
port from the National Science Foundation Graduate Re¬ 
search Fellowship under Grant No. DGE-1122492. 

Note added. —Recently, the authors of an interesting 
related experiment contacted us, alerting us to their liq¬ 
uid state NMR search for temporal order of periodically 
driven spins in star-shaped clusters [25] ■ They study a 
unique spin Hamiltonian, and they explore a range of 
cluster sizes (with N = 1,4,10, and 37 spins). 
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